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Event-by-event fluctuations in the multiplicities of charged particles and photons at 
SPS energies are discussed. Fluctuations are studied by controlling the centrality of the 
reaction and rapidity acceptance of the detectors. Results are also presented on the event- 
by-event study of correlations between the multiplicity of charged particles and photons 
to search for DCC-like signals. 



1. Introduction 

At high temperatures or high baryon number density. Quantum Chromo-Dynamics 
(QCD) describes a world of weakly interacting quarks and gluons very different from the 
hadronic world in which we live. This suggests the possibility of a phase transition as 
the temperature or density is increased in which there is a transition from a state of 
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matter where quarks are confined inside hadrons to one where quarks are free to move 
around {deconfined ) within a large volume - the quark gluon plasma (QGP). This can 
be addressed through experimental studies involving relativistic heavy-ion collisions. Ex- 
perimental searches have focused on isolating signatures of two types of phase transitions 
which might occur in extremely hot and/or dense nuclear matter. One is related to the 
deconfinement of quarks while the other is related to chiral symmetry restoration. There 
is also an interesting possibility of the existence of a tri-critical point in the phase dia- 
gram [|l|, where the transition changes from first to second order. At the tri-critical point, 
one would observe singularities in several thermodynamical variables, such as the specific 
heat and matter compressibility. The thermodynamical variables are related to event-by- 
event fiuctuations in experimental observables like particle multiplicity, transverse energy, 
and mean transverse momentum. For example, the total heat capacity is related to trans- 
verse momentum fiuctuations and matter compressibility to multiphcity fiuctuations 0,0]. 
Hence the existence of the tri-critical point can be probed in an experiment by varying the 
control parameters of the reaction (impact parameter, rapidity acceptance, beam energy, 
and system size) and studying the event-by-event fiuctuations in the global observables. 
It is also believed that in high energy heavy-ion collisions there is a possibility of creating 
a chiral symmetry restored phase. One of the possible interesting consequences of chiral 
symmetry restoration, is the formation of disoriented chiral condensates (DCC) 0. The 
detection and study of a DCC state is expected to provide valuable information about 
the chiral phase transition and vacuum structure of strong interactions. The probability 
distribution of the neutral pion fraction in a DCC domain follows the relation : 

P(/) = 1/2/7 where f = Ko/N^. (1) 

Since the majority of charged particles consist of charged pions and majority of photons 
originate from 7r° decays, DCC formation in a given domain would be associated with 
large correlated event-by-event fiuctuations in the multiplicities of charged particles and 
photons. Here we present the experimental results of multiplicity fiuctuations and correla- 
tions in relativistic heavy-ion (Pb-|-Pb) collisions as measured by the WA98 experiment 
at the CERN SPS. 



2. Multiplicity fluctuations 

Recently, much theoretical interest has been generated on the subject of event-by-event 
fiuctuations, primarily motivated by the nearly perfect Gaussian distributions of several 
observables (A^^, Nch, Et, and < pr >) for a fixed centrality bin One may define 

the relative fiuctuation as, 

where aj^ is the variance of the distribution. The fiuctuation is then studied by varying 
the control parameters of the reaction, such as the centrality and rapidity acceptance. 

The WA98 experiment at the CERN SPS has carried out a detailed analysis of photon 
and charged particle multiplicity fiuctuations. The photon multiplicity is measured 
using the photon multiplicity detector (PMD) which has the pseudo-rapidity cover- 
age from 2.9 to 4.2. The charged particle multiplicity is obtained from the silicon pad 
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multiplicity detector (SPMD) p with pseudo-rapidity coverage from 2.35 to 3.75. The 
centrality of the reaction is defined through the measurement of transverse energy (Et) 
with the mid-rapidity calorimeter (MIRAC) having a pseudo-rapidity coverage from 3.5 
to 5.5. 

2.1. Centrality dependence of multiplicity fluctuations 

It is very important to control the centrality selection carefully in fluctuation studies so 
that the impact parameter fluctuations are kept to a minimum and the distributions are 
good Gaussians. With this in mind we have used narrow centrality selections with bin 
widths of 2% in cross section, as determined from the measured total transverse energy, 
such as, 0-2%, 2-4%, 4-6%, 6-8%, etc. The fluctuations calculated for these centrality bins 
using Eqn. 1 are plotted in Fig. |I] as a function of the number of participants. The results 
from data for both charged particles and photons are compared to those obtained from the 
VENUS event generator and a simple participant model 0. In the participant model, 
where the total particle multiplicity in an event is the sum of the number of particles 
produced by each participant, the multiplicity fluctuations, cjat, can be expressed as 

UN = UJn+ {n)uJNp,rt (3) 

The value of ujn, the fluctuation in the number of particles falling within the detector 
acceptance (n) produced per participant can be obtained from nucleon-nucleon data. 
The value of ujNpart^ the fluctuation in the number of participants, can be obtained from 
simulations. From Fig. |I| we observe that the measured fluctuations are reasonably well 
reproduced by statistical models. 

2.2. Acceptance dependence of multiplicity fluctuations 

Fluctuations in the particle multiplicity have also been studied by varying the rapidity 
acceptance. The results are shown in Fig. |^. One observes, contrary to naive expec- 
tation that the fluctuations should increase as the multiplicity decreases with decrease 
in acceptance, that instead the fluctuations decrease. This can be explained through a 
simple statistical model, based on the assumption that particles are accepted following a 
binomial distribution. The details of this model are given in Ref As per this model, 
the fluctuations ) in a small acceptance region are related to the fluctuations (cUm) in 
a larger acceptance as 

UJn = l- f + fuJm (4) 

where / is the ratio of the average number of particles in the smaller acceptance to the 
average number of particles in the larger acceptance. The results from the model shown 
in Fig. 1^ are found to be in reasonably good agreement with the data. 

3. A^^ vs. Nch fluctuations 

The main motivation to study event-by-event fluctuations in the photon vs. charged 
particle multiplicity is to search for possible formation of DCCs. For this one would like 
a photon multiplicity detector and a charged particle multiplicity detector with as large a 
common coverage in 77 — phase space as possible. The main DCC search result from the 
SPS has come from the WA98 experiment |llO|JTl[| . First we discuss the results for top 5% 
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Figure 1. Multiplicity fluctuations of charged particles and photons for various centrali- 
ties in a region of common coverage of the PMD and SPMD as a function of the number 
of participants. The 7 — like results correspond to fluctuations of the measured photons 
without correction for the photon counting efficiency or purity of the photon sample. 
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Figure 2. Multiphcity fluctuations of photons and charged particles for two rj acceptance 
selections. The open squares represent estimated fluctuation values in 0.5 unit of 5r] from 
the observed fluctuations in 1.0 unit of drj. 
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Table 1 

Types of mixed events and the types of fluctuations probed. 





PMD 


SPMD 


Type of fluctuation probed: 


Ml 


Mix hits 


Mix hits 


Correlated + Individual 


M2 


Unaltered hits 


Unaltered hits 


Correlated 


M3-7 


Unaltered hits 


Mix hits 


A^-^ only 


M3-cli 


Mix hits 


Unaltered hits 


Nch only 



most central events. For the pseudo-rapidity and azimuthal region common to both PMD 
and SPMD used for this analysis, the average photon and charged particle multiplicity is 
335 and 323, respectively. The experimental results are compared to simulated events and 
various types of mixed events to investigate possible DCC formation. The simulated events 
were generated by passing the VENUS output through a detector simulation package 
(GEANT), which incorporates the WA98 experimental setup. Four different types of 
mixed events were generated from the data with equivalent event sample sizes as for real 



events. The details of the construction of the mixed events can be found in Ref. [|T0|. 
In Table 1 we summarize the construction and the physics issue probed by each type 
of mixed event. A common analysis was carried out on the data, simulated, and mixed 
events to investigate the source of any observed fluctuations. 

3.1. A^^ vs. Nch correlation analysis 

Analysis of the correlation between and A^^ch is useful to search for DCC-type fluctu- 
ations. The details of this can be found in Ref. P,pi]|. From the event-by-event correlation 
between A^^-uke and A'ch in various 0-segments (obtained by dividing the 0-space into 2, 
4, 8, and 16 bins) and a common correlation axis (Z), one can obtain the closest distance 
(Dz) of the data points to the correlation axis. The correlation plot is shown in Fig. ^. 
In order to compare the fluctuations for different bins on a similar footing, a scaled 
variable, Sz = Dz/s{Dz), is used, where s{Dz) represents the RMS deviation of the Dz 
distribution for VENUS events analyzed in the same manner. The width of the distri- 
bution of Sz represents the relative fluctuations of A'^-ukc and A^ch from the correlation 
axis at any given (f) bin. Since the width of the Sz distribution quantifies the amount of 
fluctuation, the RMS deviations of these distributions for data are compared with those 
from mixed events and simulations in Fig. |. 

The RMS deviations of the Sz distributions for mixed events and data agree for 1 bin 
in (by construction) and for 16 bins in within the quoted errors. However, for the 
other bins in one observes that results from Ml mixed events are lower than those from 
data. This indicates the presence of localised non-statistical fluctuations in the data. The 
source of the additional fluctuations is understood by the comparison to the M3-type of 
mixed events. Comparison shows that the excess fluctuations have contributions both 
from A^^ and N^h- However the RMS deviations of the M2-type of mixed events closely 
follow those from data, suggesting the absence of event-by-event correlated DCC-like A^^ 
vs. Nch fluctuations. 
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Figure 3. The event-by-event correlation between Net and A^^_iikc for the top 5% cen- 
trahty class. Overlaid on the plot is the common correlation axis (Z-axis). 
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Figure 4. The root mean square (RMS) deviations of the Sz distribution for various 
divisions in the azimuthal angle. Errors include both statistical and systematic sources. 
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Figure 5. The root mean square (RMS) deviations of the FFC distribution for various 
divisions in the azimuthal angle. Errors are due to both statistical and systematic sources. 



3.2. Multi-resolution analysis based on discrete wavelet transformation 

A multi-resolution analysis using discrete wavelet transformations (DWT) has been 
shown to be quite powerful in the search for localized domains of DCC |]12[. For the 
present DWT analysis the full azimuthal region is divided into smaller bins in 0, the 
number of bins at a given scale j being The input to the analysis is a spectrum 
of the sample function at the smallest bin in corresponding to the highest resolution 
scale, imax{= 5 here). In the present case, the sample function is chosen to be the photon 
fraction, given as, 

r(0) = iV^-iikc(0)/(iV7-iikc(0) + iVch(0)) (5) 

The output of the DWT consists of a set of wavelet or father function coefficients (FFC) 
at each scale, from j = l,...,{jmax — !)• The coefficients obtained at a given scale, j, are 
derived from the distribution of the sample function at one higher scale, j + 1. The FFCs 
quantify the deviation of the bin-to-bin fluctuations in the sample function at that higher 
scale relative to the average behavior. The presence of localized non-statistical fluctuations 
will increase the RMS deviation of the distribution of FFCs and may result in non- 



Gaussian tails [T^ . Once again, comparing the RMS deviations of the FFC distributions 
of data, mixed events, and VENUS events may allow to draw inference about the presence 
of localized fluctuations. 

The RMS deviations of these FFC distributions are summarized in Fig. ^. The RMS 
deviations of the FFC distributions for the data, VENUS, and mixed events are found 
to be close to each other (within quoted errors) for the case of 32 bins in 0. While 
the values for M2 mixed events are found to closely follow those of the data for all bins 
in 0, the RMS deviations for the M3 mixed events lie between those of the data and 
Ml mixed events. These results are consistent with those obtained from the analysis of 
the Sz distributions. These observations indicate the absence of event-by-event localized 
correlated fluctuations (DCC-like) between iV^-uke and A^'ch- However, they do indicate 
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the presence of localized independent fluctuations in both photon and charged particle 
multiplicities for intermediate bin sizes in azimuth. 

3.3. Centrality dependence of A^^ vs. Nch fluctuations 

It is interesting to study the centrality dependence of A^-^ vs. Nch fluctuations |jl3l . This 



study has been carried out for four centrality classes corresponding to the top 5%, 5% 
- 10%, 15% - 30%, and 45% — 55% of the minimum bias cross section. The correlation 
and DWT analysis was carried out on events from each of the centrality classes and the 
results in terms of RMS deviations of Sz and FFC distributions for data, various sets of 
mixed events, and simulations were compared. For all four classes, the RMS deviations 
of the Sz and FFC distributions of data and M2 mixed events agreed reasonably well 
with each other. However, comparison to Ml and M3 mixed events showed the presence 
of localised uncorrelated event-by-event non-statistical fluctuations in both photon and 
charged particle multiplicities. In order to quantify the strength of the total localised 
A^-y_iike and N^h fluctuations for various bins in cf) and for different centrality classes, we 
define a sensitivity parameter x as : 



X=^^ (6) 

where si and s correspond to the RMS deviations of the FFC distributions of the Ml mixed 
events and real data, respectively. The results are shown in Fig. ^ as a function of the 
number of bins in cf) for the four different centrality classes. Qualitatively similar results 
are obtained when x is calculated using the RMS deviations of the Sz distributions. 
The shaded portion indicates the region of x where s is one a greater than the RMS 
deviation FFC distributions for Ml events, where a is the total error on the Ml event 
RMS deviation. It represents the limit above which a signal is detectable. The result 
shows that the strength of the fluctuations decreases as the number of bins in increases, 
with a strength which decreases to below detectable level (within the quoted errors) for 16 
and 32 bins. It is also observed that the strength of the signal decreases with decreasing 
centrality for 4 and 8 bins in azimuthal angle, although the tendency is not very strong. 

3.4. Upper Limit on DCC production at SPS 

If the amount of DCC-like fluctuations in the experimental data were large, then the 
RMS deviations of Sz and FFC distributions for data would have been larger than those 
of M2 events. Since this is not the case, we have extracted upper limits on the probability 
of DCC-like fluctuations at the 90% confidence level following the standard procedure 



as discussed in Ref. |TT],|T^. This is done within the context of the simple DCC model 
described in Ref. ||ll|. We give the upper limits for the two most central event classes 
(0-5% and 5-10%). The 90% CL upper limit contour has been calculated as % + 1.28e-^, 
where x is calculated using Eqn.(^. Here si and s correspond to the RMS deviations of 
the FFC (or Sz) distributions for M2 mixed events and real data, respectively, and is 
the error in x from the FFC (or Sz) analysis. To relate the measured upper limit on the 
size of the fluctuations in terms of DCC domain size and frequency of occurrence we take 
use the simple simulated DCC model described in Ref. The upper limit is set at that 
value of frequency of occurrence for a fixed DCC domain size at which the x value from 
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Figure 6. The fluctuation strength parameter for the four centrahty classes. The error 
bars are shown only for the top 5% centrality class for clarity of presentation. The errors 
are similar for the other centralities. 



the DCC model matches with that of the x + l-28e^ upper limit from the experimental 
data. The results for centrality classes 1 and 2 are shown in Fig. 0. 

4. Summary 

A detailed study of centrality and rapidity acceptance dependence of multiplicity fluc- 
tuations carried out by the WA98 experiment at the CERN SPS shows an absence of 
significant non-statistical fluctuations in the photon and charged particle multiplicity. A 
model independent study of event-by-event correlated fluctuations (DCC-type) in photon 
and charged particle multiplicity using a robust mixed event technique reveals an absence 
of significant DCC- like correlated fluctuations at the SPS. However, the analysis indi- 
cates the presence of event-by-event uncorrelated fluctuations in both A^^ and A'ch beyond 
those observed in simulated and mixed events for limited region of azimuthal angle, with 
the strength of fluctuation increasing with increase in centrality. Using the results from 
the data, mixed events, and a simple model of DCC formation, an upper limit on DCC 
production in Pb-|-Pb collisions at SPS energies has been set. 
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Figure 7. The 90% confidence level upper limit on DCC production for central Pb+Pb 
collision at 158 A GeV/c as a function of the DCC domain size in azimuthal angle. The 
solid line corresponds to data from the top 5% and dashed line to top 5 — 10% of the 
minimum bias cross section as determined by selection on the measured transverse energy 
distribution. 
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